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We explore in detail oscillations of the solar ^Be neutrinos in the matter of the Earth. The depth of 
oscillations is about (0.1 — 0.2)% and the length ~ 30 km. The period of the oscillatory modulations 
in the energy scale is comparable with the width of the line determined by the temperature in the 
center of the Sun. The latter means that depending on the length of trajectory (nadir angle) one 
obtains different degree of averaging of oscillations. Exploring these oscillations it is possible to 
measure the width of the ^Be line and therefore the temperature of the Sun, determine precisely 
Amii, perform tomography of the Earth, in particular, measure the deviation of its form from 
sphere, and detect small structures. Studies of the Be neutrinos open up a possibility to test 
quantum mechanics of neutrino oscillations and search for the sterile neutrinos. Accuracy of these 
measurements with future scintillator (or scintillator uploaded) detectors of the ~ 100 kton mass 
scale is estimated. 

PACS numbers: 14.60.Pq, 26.65.+t, 95.85.Ry, 95.60.Jw, 


INTRODUCTION 

Exploration of the solar neutrinos moves to the phase 
of precision measurements. In future, developments of 
the detection techniques and construction of the large 
mass detectors can open up various new possibilities and 
one of them is study of the solar ^Be neutrino oscillations 
in the Earth. 

The solar ^Be neutrinos have two salient properties: 
low energies and narrow width. Due to environment ef¬ 
fects in the Sun [T] the energy profile of the ^Be neutrinos, 
g{E), is broadened and has an asymmetric form. The 
maximum value of g{E), E ~ 862.27 keV, is shifted to 
higher energy with respect to the laboratory value Eiat ~ 
861.64 keV [2]. The left (low energy) wing of the profile 
is determined by the Doppler shift caused by thermal ve¬ 
locities of ^Be nuclei. The right (high energy) wing is 
mainly due to collisions determined by temperature at 
the center of the Sun (averaged over ^Be neutrinos pro¬ 
duction region). The width of the line at the half of 
height equals 

Tse ^ 1.6 keV, (1) 

so that the relative size of the width is 

% = 1 . 86 - 10 "®. ( 2 ) 

h/ 

The width and the shift of maximum of g{E) are pro¬ 
portional to the central temperature of the Sun Tc- In [T] 
it was proposed to determine Tc by measuring the shift of 
the peak. It was also mentioned that finite width of the 
^Be spectrum affects the depth of vacuum oscillations. 

According to the standard solar model the flux of ^Be 
neutrinos is known with 1.4% accuracy. For the LMA 


MSW values of oscillation parameters the ^Be neutrino 
line is at the low energy edge of the so called transi¬ 
tion region between the averaged vacuum oscillations and 
matter dominated conversion. The flux is suppressed 
mainly by the averaged vacuum oscillations probability 
with small additional suppression due to matter effect in 
the Sun. 

The ^Be neutrino flux measured by BOREXINO [3] , [1] 
and KamLAND [5] is in a good agreement with the SSM 
predicted flux with suppression according the LMA solu¬ 
tion. In fact, these measurements (reached 5% accuracy) 
provided important confirmation and consistency checks 
of the LMA MSW solution of the solar neutrino problem. 
Time dependence of the flux has been explored [5]. The 
data are in agreement with seasonal variations due to ec¬ 
centricity of the Earth orbit and no other variations have 
been found in agreement with LMA MSW expectations. 
In particular, after 2 years of exposure and 0.133 kton 
fiducial mass the following bound on the Day - Night 
asymmetry has been obtained [^: 


Adn = 2 ^—= 0.001 ± 0.012(stat) ± 0.007(syst). 


On the way from the Sun the coherence of neutrino 
state is lost, so that incoherent fluxes of the mass states 
arrive at the surface of the Earth. Inside the Earth these 
mass states oscillate due to matter effect. According to 
LMA MSW solution these oscillations proceed in the low 
matter density regime, and the expected effects are very 
small (see e.g. [7], [B] [3], [IH], [U, HU, |I3j, [H])- The 



effects are determined by the parameter: 


2. OSCILLATION IN THE SUN AND THE 
EARTH 
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where E is the neutrino energy, 14 = ^/^Gpne is the mat¬ 
ter potential with Gp and rig being the Fermi coupling 
constant and the electron number density correspond¬ 
ingly. The parameter e characterizes deviations of the 
mixing angle and the oscillation length in matter from 
their vacuum values: 9i2 and 


L = 


AttE 

Am|i 


28.5km 


7.5 X IQ-^eV^ 

Am|i 


( 5 ) 


In fact, e determines the depth of oscillations of neutri¬ 
nos with definite mass, « 0.5e which is about ^ 0.1% 
in the mantle and ~ 0.2% in the core. Being of the or¬ 
der 10 “^ the expected effects are far beyond the present 
BOREXINO ^ as well as expected SNO+ sensitivities. 

Next generation of large (several tenth of ktons to 
hundred ktons) scintillator detectors like JUNO [IS] or 
LENA [16] will have sub-percent sensitivity to the Day- 
Night asymmetry. Higher sensitivity can be achieved 
with 100 kton mass scale scintillator uploaded water de¬ 
tectors, WBLS [T7]. For 100 kton fiducial mass and 5 
years exposure such a detector will collect 1.9-10^ bigger 
statistics than the one used for the result Corre¬ 
spondingly, the statistical error will be reduced down to 
3 • 10“^. So, if systematic errors is well controlled, the 
0.1% size Earth matter effects on the ^Be neutrinos can 
be established at about 3cr level. 

In this connection we will explore in detail the ^Be neu¬ 
trino oscillations in the matter of the Earth. There are 
two very interesting coincidences related to the energy 
and width of the ^Be neutrino profile which allow one 
to obtain in principle unique information about neutrino 
properties, characteristics of the ^Be neutrino spectrum, 
properties of the Earth density profile and quantum me¬ 
chanics of neutrino oscillations. We estimate possibili¬ 
ties of future large detectors to determine the width of 
the line, and consequently, the central temperature of 
the Sun), to measure Am^i with unprecedent accuracy, 
to perform tomography of the Earth and search for very 
light sterile neutrinos. 

The paper is organized as follows. In Sec. 2 we present 
relevant analytic results for the probabilities of oscilla¬ 
tions in the Sun and the Earth as well as compute the 
relative variations of the flux with the nadir angle. In 
Sec. 3 the effects of averaging of the flux over the ^Be 
energy spectrum are explored. We compute the time 
(nadir angle) variation of number of events and estimate 
a potential of future 100 kton scale detectors to establish 
the Earth matter effect, to measure Fee and in 

Sec. 4. Searches for sterile neutrinos are considered in 
Sec. 5. We conclude in Sec. 6 . 


The electron neutrino Ve produced in the center of the 
Sun is adiabatically converted into the combination of 
the mass eigenstates vi, V 2 , The combination is deter¬ 
mined by the mixing angles, 0^2 = 
in the production point 

Pe Ci3 cos 9 ^ h'l -I-Ci3 sin05^ 1/2 + Slat's, (6) 


where C 13 = cos 6 * 13 , S 13 = sin 6 * 13 , and we neglected the 
matter effect on 1-3 mixing, so that 9^{V^) «013 m- 
We use the standard parametrization of the PMNS mix¬ 
ing matrix. In Eq. 0)^ is the value of 1-2 angle in 
matter averaged over the ^Be neutrino production region. 
It is given by 

tan2 20™ « tan^ 20i2 (1 - , (7) 

V Am|i cos 2012 / 

where is the average matter potential in the produc¬ 
tion region. 

On the way from the Sun to the Earth the wave pack¬ 
ets of mass eigenstates Pi spread and separate. At the 
production the lengths of the wave packets in the config¬ 
uration space equal 

Stt 

~ 4 ^— = 6 • 10 “® cm. ( 8 ) 

i Be 


Due to spread their sizes become at the surface of the 
Earth 


^spread 


^su'nX'Be 


(9) 


where is the absolute value of mass of Vi, and Lsun is 
the distance from the Sun to the Earth. For hierarchical 
spectrum, m 2 « •\/Am|)', we obtain from (j^ ^spread _ 
3.8 • 10“® cm, (j.e. 2 orders of magnitude larger than the 
original size of the packet). However, separation of the 
wave packets of different mass eigenstates is larger: 


Aa; = Ls„„Am 2 i/ 2 A^ = 2-10 ^ cm. 


Even in the case of degenerate spectrum, mi « m 2 « 0.1 
eV, we obtain « 0.25. In addition there is 

also averaging of oscillations over the neutrino production 
region inside the Sun. 

Due to loss of coherence, neutrinos arrive at the sur¬ 
face of the Earth as incoherent fluxes of vi, V 2 and 1/3 
with relative admixtures given according to Eq. Q by 
Ci 3 cos^ 9^2, Ci 3 sin^ 0)^ and S 13 correspondingly. In mat¬ 
ter of the Earth each of these mass states splits into eigen¬ 
states in matter and oscillates. 

The probability to find t'e in the detector after crossing 
the Earth can be written as 


P = C?3(cOs2 0™Pie + sin2 0™P2e) + 4 

= c? 3 (cos 20(^Pie -f c ?3 sin^ 0^^) -f sfs, (10) 
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where Pie and P 2 e are the probabilities of vi —> and 
1^2 Ve transitions in the Earth correspondingly. Here 
we used the unitarity relation: Pie + P 2 e + S 13 = 1 or 
Pie + P2e = 0 ^ 3 . 

The probability P can be represented as 

P = Pd + AP, ( 11 ) 

where 

(A 

P£) = A^(l + cos 20)^ cos 26*12) + Sj3 (12) 

is the probability during the day when Pie = Pfe = 
C?3 COS^ 6*12, P2e = P2°e = <^3 6*12, and 

AP = els cos 20™ (Ae - P°e) (13) 

is the difference of the probabilities during the day and 
the night. The probability Pie equals [5] 


Combining Eqs. (12), (13) and (14) we can write the 
relative variation of the Vg—Oxiyi due to the Earth matter 
effect as 


^0 = 


AP 


=-^/(Ato^i,012,0i3) J dxVe{x)sin(j)^_^j^, 

(19) 


where 


/(Amli, 0 i 2 ) = 


2 cos 201 ^ sin^ 20 i 2 


1 + cos 20 )^ cos 2012 + 2 tarn 0 i 3 


( 20 ) 

For 012 = 34° and Ami 2 = 7.5 x 10 ® eV^ we obtain 
/(Amii, 0 i 2 ) ~ 0 .43; corrections due to the 1-3 mixing 
are below 0 . 1 %. 

Quick estimation of the effect can be done for con¬ 
stant density profile. In this case we obtain = 

A()j(L — x), and in the lowest order in e the Eq. (19) 
gives 


Ple/C 


2 _ 
13 — 


where 



dx Ve{x) sin 

(14) 


= ^dx Arnix), (15) 


/\ yT} ^ j 

Am{x)= Y [cos 2012 - clse{x)Y + sin^ 20 i 2 (16) 

and e is determined in Here x is the distance be¬ 
tween an entrance point to the Earth and a given point 
of trajectory. 

The oscillation length can be written as 

O-TT 

= w— -j-T = [1 + COS20i2Ci3e -|- 0(e^)]. (17) 


The matter correction, cos 20 i 2 Ci 3 e, is of the order (1 — 
2 ) • 10-3. 


Notice that since in the Sun P > 1/2, or cos 20)^ > 0, 
the oscillations in the Earth suppress the survival prob¬ 
ability, AP < 0, in contrast to high energy Boron neu¬ 
trinos for which cos 20 )^ < 0 and partial regeneration of 
the Vf. flux occurs. 

The expression in (14) is equivalent to the result of 
adiabatic perturbation theory m- 


Pie/cls = cos^ 012 - sin20*2 sin2(032 - 012) Sin^ (/)o-!-l /2 

-H sin2(2032 - 012)/ de^ix) cosc/™.,.^, (18) 

Jo 


where 0 f 2 is the mixing angle in matter at the surface of 
the Earth. The angle 0 f 2 is given by Eq. ([^ with V/ sub¬ 
stituted by Veipsurf)- We use this formula for numerical 
computations inside different layers of the Earth. 


Ki-clsef siii^ ^AmL. ( 21 ) 

So, the depth of the oscillations equals chef « ef. 

Dependence of AP/P” on the nadir angle p is shown 
in Figs. (red lines). In our computations we used 
Pie from Eq. ( [T8| , the spherically symmetric Earth and 
the 5 layers parametrization of the density profile [5T] 
with sharp density jumps at 410 km, 660 km, 2830 km 
and 5150 km from the surface of the Earth. The length 
of the neutrino trajectory inside the Earth is given by 
L = 2ReCOSp; so dependence on cosp is equivalent to 
the dependence on L and therefore Figs. and reflect 
dependence of oscillations on distance. Trajectories with 
cos 77 > 0.83 cross the core of the Earth. 

According to the Figs. 0 |(red lines) the depth of os¬ 
cillations changes with the nadir angle which is related to 
breaking of adiabaticity. If the adiabaticity condition is 
fulfilled along whole the trajectory, the oscillation depth 
would be determined exclusively by the mixing angle at 
the surface of the Earth. (More precisely, by the mixing 
angle averaged over the distance of the order of oscilla¬ 
tion length). This is satisfied for shallow trajectories with 
77 >1.21. At 77 « 1.21 neutrinos cross the first density 
jump (410 km from the surface) and at 77 « 1.11 - the 
second one (660 km), where adiabaticity is broken, and 
consequently, the depth of oscillations increases. 

The biggest change is at 77 = 0.58, where neutrinos 
start to cross the core of the Earth. The depth of oscilla¬ 
tions increases by factor of 2. At 77 = 0.19 neutrino tra¬ 
jectories cross also borders of the inner core. Breakdown 
of the adiabaticity and oscillations in different layers lead 
also to modulations of the oscillatory picture which is re¬ 
lated to the interference of effects from different layers 
(see below). 

The length of the deepest trajectory in the mantle is 
Lmax = 10570 km, therefore the number of periods in 
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the mantle range, rj > 0.58, equals 

« 370. (22) 

t'l/ 

The number of periods for the core crossing trajectories 
is bigger: about 420 for the vertical direction. 


3. EFFECT OF THE BE NEUTRINO LINE 
WIDTH 

Averaging over the energy 

Signal in the detector is determined by the variation, 
Ag, averaged over the energy spectrum of the ^Be neutri¬ 
nos. Energy resolution of a detector is much worse and 
therefore irrelevant for the oscillation picture. The ^Be 
neutrino flux can be written as 

FBe{E,t) = FUt) 9(E), (23) 

where Eg^{t) is the total flux and g{E) is the energy 
profile normalized to 1: J dEg{E) = 1 . Dependence of 
Fgg on time is due to the eccentricity of the Earth orbit 
around the Sun, as well as due to rotation of the Earth 
itself (due to bigger distance from the Sun at night). The 
former (±3%) must be taken into account, the latter (< 
10 “'^%) can be neglected. 

Variation of the z/g—flux averaged over the ^Be spec¬ 
trum for the spherically symmetric Earth equals 

A,{t) = J dE g{E) A,{E,t). (24) 

The dependence of Af. on t] is shown in Eigs. (blue 
solid lines). In what follows we will present complete 
interpretation of the this dependence. 

1. Suppression of the oscillation depth. Integration 
over the energy profile leads to partial suppression of the 
oscillation depth which depends on the nadir angle. At 
the same time the integration does not affect the average 
value of the relative variations. The change of the oscil¬ 
lation depth with the nadir angle is result of interesting 
and accidental coincidence: the width of the line, Tne, 
is of the order of the period of oscillatory curve in the 
energy scale, /S.Et'. 


AEt ~ TBe- 


(25) 


which depends on rj. Therefore relation between AEt 
and Tee, and consequently, change of the degree of av¬ 
eraging and the depth of oscillations depend on rj. For 
shallow trajectories (large rj) AEt ^ Tee the averaging 
effect is negligible. With decrease of p the period AEt 
decreases and averaging becomes stronger (see Fig. [^. 
For the deepest trajectories in the mantle the suppression 
is characterized by factor 2-3. 

In configuration space, partial averaging corresponds 
to partial lost of coherence during propagation inside the 
Earth. Indeed, the relative shift of the wave packets of 
the eigenstates equals 


Ax = L—^ = 5-10“®cm 


(—) 

\ lO^km J 


= O.ScTo 




\ I0"^km J 
(28) 

In the last equality in (28) we used estimation (|^. The 
shift is of the order of original size of the wave packet. 
Ax ~ ax, which is complementary to the relation (25). 


Since separation of wave packets is partial the effect of 
averaging depends on the shape of the wave packet. So, 
measuring dependence of the depth of variations of signal 
on the base-line ( 77 , see Fig. 2) one can in principle restore 
the shape of the wave packet. (Loss of coherence is due 
to relative shift of the wave packets and the coherence 
condition is determined by the original size of the wave 
packet and not by the length of the packet after spread 
m)- This is very rare situation; another one is realized 
for supernova neutrinos propagating in the Earth m- 


2 . Enhanced effect for the core crossing trajectories. 
For trajectories with 77 < 0.58 the depth of oscillations 
averaged over the neutrino energy is factor of 4 - 5 larger 
than for the deepest trajectories in the mantle. Further¬ 
more, AP is positive during some part of the oscillation 
period. Both effects are due to higher matter density 
in the core and interplay of oscillation effect in the core 
and mantle. Appearance of the positive Earth matter 
effect, AP > 0, is also related to that fact that the av¬ 
erage values of A^ in the core and in the mantle should 
be the same. The latter can be explained by the atten¬ 
uation effect (see below): with strong averaging the core 
should not be seen. The only what one can observe is 
the average probability which is the same in the core and 
mantle. 


The period is defined by equality 

AEr = 27r, 


dE 


(26) 


where 4)^ ~ ^irL/lm ~ 27rL/l^ is the oscillation phase. 
We find from 


AEt =E— =E 


L 


L 2RECOsr]’ 


(27) 


3. Attenuation effect. Integration over the energy leads 
to lost of sensitivity to the remote structures of the den¬ 
sity profile [5]. This is characterized by the attenuation 
length, Xatt, in such a way that structures at distances 
from a detector bigger than Xatt are not resolved in the 
oscillation picture. To find the dependence of Xatt on 77 
explicitly we take for g{E) a Gaussian profile with the 
width TBe- (The use of exact line shape does not change 


this result substantially.) In this case the expression (24) 
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FIG. 1: The relative variations of the electron neutrino flux as function of the nadir angle of the neutrino trajectory. Dotted 
(red) line shows without averaging; solid (blue) line is Ae which corresponds to the variations averaged over the energy 
spectrum of the ^Be neutrinos. Spherical symmetry of the Earth is assumed. 




FIG. 2: The same as in Fig. [^with two zoomed regions: r] = [0.58 — 65] (left) and [1.3 — 1.35] (right). 


gives 


Ae =-^c?3/(A?7^^6») J dx V{x)F{L-x) sin 

(29) 


where 


F{L — x) = exp < —2 


7rrBe(i - x) 


1 2 


(30) 


According to Eqs. (29) and (30) a contribution to the 
integral (29) from structures situated at the distance {L— 
x) > Xatt from detector, where the attenuation length 
Xatt is determined by the relation 

^ TrTBeXatt ^ 

EL 


becomes exponentially suppressed. This gives Xatt 
7000 km which is comparable with radius of the Earth 
and bigger than the distance to the core. So, the core of 


the Earth should be seen in agreement with result of Fig. 
[^(increase of the oscillation amplitude at rj = 0.58). The 
average values of oscillation curve in the mantle and the 
core are nearly the same. 

4. Modulation of the oscillation picture; interference. 
At r] = 1.21 and ry = 1.11 the oscillatory picture changes: 
the depths of oscillations slightly increases, and high fre¬ 
quency oscillations start to be modulated by low fre¬ 
quency modes. Periods of modulations increase with de¬ 
crease of rj (see Fig. [^. The modulations originate from 
the density jumps at distances hi = 410 km (by 5%) and 
at ^2 = 660 km (by 10%) from the surface of the Earth. 
These density changes occur over depth interval less than 
5 km m, m, which is much smaller than the oscilla¬ 
tion length. Therefore the jumps break adiabaticity. The 
adiabaticity is well satisfied within the layers bounded by 
the jumps. The jumps lead to increase of the depth of 
oscillations at sinry = 1 — hi/R e- The periods of mod¬ 
ulations in the ry scale can be obtained in the following 
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way. The length of the trajectory with a given value of 
the nadir angle r] between the surface of the Earth and 
the ith jump equals 

Xi{r]) = Re COST] — R\ cos^ t] — 2REhi + (31) 


Consequently, the length of trajectory in the layer be¬ 
tween the i-th and j-th jumps is Xij = Xi^rj) — Xj(ri). 
Then the period of modulations, 6{r]) = rji — r ]2 is deter¬ 
mined by the condition 


[xiim) - Xjim)] - [xiim) - Xj{ri2)] = im- (32) 


In the approximation of flat layers {Re —t oo) or 
Re cost] ^ hi the condition simplifies: the length of the 
trajectory equals approximately hij cos 77 and instead of 
(321 we obtain 


A(cos 77 ) = (cos 771 — cos 772 ) « ^ cos^ 77 . (33) 

tli 


From ( [^ [M| ) for 77 = 0.58 (the deepest mantle trajec¬ 
tory) we find A 77 = 0.08 in agreement with the results 
of Fig. With increase of 77 the period of modulations 
decreases. Thus, characteristics of modulations encode 
information about parameters of non-adiabatic density 
jumps (their position and size). Modulations become es¬ 
pecially profound for the deepest mantle trajectories. 


Geophysics of the Be neutrino oscillations 


At the sub-percent level of the experimental accuracy 
and very short oscillation length a number of new effects 
become accessible and should be taken into account. Re¬ 
call that the oscillatory curves shown in Figs. [Hi corre¬ 
spond to the ideal spherically symmetric Earth. In reality 
the Earth profile is not symmetric and deviations are of 
the order of oscillation length. Deviations include: 

(i) Non-sphericity of the Earth as whole. 

(ii) Small scale structures at the surface (mountains, 
oceans, etc.). 

(iii) Structures in the crust (layers with anomalous 
density, cavities). 

(iv) Possible structures in the mantle and the core, in 
particular deviation of the core from spherical form. 

Effects of the deviations on can be quantified by 
contributions to the total oscillation phase which can 
be written according to (15), (16) as 


$ 


m 


1 — cos 20 i 2 Ci 3 ^ J dxe{x) 


27r^(l — cos20i2Cj3e). 


(34) 


Here e(L) is the averaged value of the matter parame¬ 
ter over the neutrino trajectory with the length L. The 


matter correction to (last term in (34)) equals 
L 


A4.„ = 2,-co.2e„f=|(j^) 


( — ) 

I 7000km J 


. (35) 


For 77 = 1.2 (the outer mantle) we obtain A4>m = O.Ott. 
For 77 = 0.58 (the deepest trajectory in the mantle) the 
corresponding number is A$m = Stt. For 77 = 0.4 (the 
core crossing trajectory) we find Ad)^ = 5 . 97 r. So, to 
avoid an uncertainty in the phase, e should be known 
with accuracy better than 20%, 8 % and 4% in the outer 
mantle, deep mantle and core crossing trajectories. In 
general, all possible corrections and uncertainties (e.g. 
in L, Am 2 i, etc.) which change the phase by ^ tt/O 
should be taken into account. In view of many oscillation 
periods obtained over the baseline even small effects may 
become important. Large uncertainties in the density 
profile would wash out the oscillatory dependence. 

Let us consider some effects of density profile pertur¬ 
bations. 

1. Deviation of the Earth shape from ideal sphere can 
be characterized as follows. The polar and equatorial 
radii equal Rpoi = 6356.7 km and R^q = 6378.1 km cor¬ 
respondingly. The difference in diameters is about 43 
km. For a given latitude tp the distance from the surface 
to the center of the Earth can be approximated by 


R^ « Reg - {Req “ Rpol)' sin^ Ip = R^q - 21.5km • sin^ ip. 

Introducing the average radius of the Earth as R = 
0.5(i?eg + Rpoi) we obtain the deviation from the average 
R^ — R = {Req — Rpoz)(0.5 — sin^ Ip). E.g. for equator 
{ip = 0 ) the difference of trajectory lengths can be as 
large as {Req — Rpoi)/2 = 21.4 km. 

Due to non-sphericity of the Earth for the same value 
of the nadir angle the length of trajectory and conse¬ 
quently the oscillation phase depend on the azimuthal 
angle. Therefore deviation from sphericity leads to mod¬ 
ification of the oscillatory curves in Figs. m If variable 
77 is used, one needs to average phase over the azimuthal 
angle. This however, may lead to complete averaging of 
the oscillations. Indeed, for a detector at the latitude ip 
and for small nadir angle 77 the difference of lengths of 
trajectories for different azimuthal angles can be as large 
as 


6x ^ 2R{r — l)tanT]sm2ip (36) 

It can be about (10 - 20) km even for small 77 . 

Inversely, due to non-sphericity the trajectories with 
the same L have different 77 and azimuthal angles and 
therefore L does not fix the density profile uniquely. So e 
can be different for trajectories with the same L but dif¬ 
ferent 77 . This variation, however, can be neglected in the 
first approximation, and one can use L as the parameter 
to mark events. 

2. The effect of small structures at the surface of the 
Earth depends on 77 . Let us neglect here non-sphericity 
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and use r] variable. Small structures produce distortion of 
the periodic sinusoidal curves shown in Figs. [pleading 
to appearance of irregular perturbations, substructures, 
shifts of maxima and minima, etc.. 

Consider perturbation produced by a structure at the 
point where neutrino enters the Earth (the beginning of 
the trajectory) which has the surface length Xpert and 
height hpert- Modification of Ae(ri) produced by this 
structure depends on relative values of Xpert and the dis¬ 
tance xt which corresponds to rjx - period of Ae(r]) in 
rj. In turn, r]T is determined from the equation 


2Re [cos rj - cos{r] + t^t)] = Im ~ lu (37) 


which can be rewritten as 

sin ?7 sin Tyr + cos 77(1 — cos 77 t) = (38) 

2Re 

Then the distance xt as function of 77 is given by 

XT ~ \/ll + [4i?B cos?7sin(0.5?7T)]^. (39) 


For very shallow trajectories (77 ~ 1.5), xt ~ ^ 1 / ~ 28 
km, which is comparable with the length of perturba¬ 
tions. 

With decrease of 77 the xt increases. For directions not 
very close to vertical the first term in (38) dominates and 
we obtain sinTyT’ « 1^/i‘2RE sin rj), or 


Vt 


^ra 

2Re sin 77 


Insertion of this r]T into (39) gives 

L 

XT 


smTy 


(40) 


(41) 


E.g. for the deepest trajectories in the mantle, sin 77 = 
0.56, we obtain xt = 50 km. xt increases faster near 


vertical directions. For 77 = 0 we find from (38 1 and (39) 


XT ~ 2\/I^Re « 850 km. 


The size of perturbation of Ae is given by the addi¬ 
tional phase acquired due to perturbation of the profile: 

^ o ^pert 

^pert ~ j 5 

where Lpert ~ hpert for nearly vertical trajectories, and 
Lpert ~ Xpert for nearly horizontal trajectories. The 
change of Ae equals 

AAe « [cos - C0S($„ -k $pert)]- (42) 

For shallow trajectories already lpert ~ Xpert = 5 km 
become important, ^pert = 0 . 357 r and it will change Ae 
in the interval A 77 ~ Xpert Ixt ~ Xpert/lu = 0.17 of the 
period. The changes of Ae in units = —cf^ef are 


0 —)■ 0.2, 0.5 —)■ 0.95, 1 —> 0.73 for = 0, 7 r/ 2 , tt cor¬ 
respondingly. This also means that maximum and min¬ 
imum will be shifted. For Xpert = 15 km the effect will 
be of the order 1 , e.g. maximum will become minimum 
and vice versa. In this case the change can be viewed as 
a local shift of the oscillatory curve by half of period in 
77 . 

With decrease of 77 effect of perturbation of the same 
size will decrease since Xt becomes larger and lpert be¬ 
comes smaller. That is, the change will be over smaller 
part of the period and the change of Ag will be smaller. 

For vertical trajectories we have Xpert ^ xt and 
^pert < 0 . 3577 . So, the effect shows up as small pertur¬ 
bations of oscillatory curve with typical size being much 
smaller than the period. 

Since the oscillation length is about 30 km, small struc¬ 
tures at the surface of the Earth (mountains, oceans, 
seas) as well as in the crust may become important. Rel¬ 
atively small structures can be averaged out when inte¬ 
grating over the nadir angle intervals. 

3. It is more difficult to control underground struc¬ 
tures. Although variations of density in regions of (5 
- 30) km size do not produce significant change of the 
phase and their effect will be subleading . In principle 
one can study the shape of the core of the Earth and its 
deviation from sphere. 

In this connection let us estimate effect of finite size of 
the ^Be neutrino production in the Sun. From the Earth 
the ^Be neutrino production region in the Sun is seen as 
a disk of the angular size 'qdisc ~ 10“^ radians. Neutrinos 
produced in different parts of the disc will have slightly 
different nadir angles and therefore baselines. Therefore 
one needs to perform integration over the corresponding 
intervals of the nadir angle. The value of r^disc is an 
order of magnitude smaller than the oscillatory period 
A 77 « 0.002 — 0.003, so the effect of averaging over disc 
is expected to be small. 

To quantify this we approximate the distribution of the 
neutrino sources in the Sun by 

= (43) 

Ufj \ 2tlT 

with Sp — 1-9 X 10“"^. The result shown in Fig.j^ indeed, 
confirms the smallness of the effect. Averaging over the 
disc is slightly stronger for large nadir angles where the 
oscillatory period is smaller. 


4. OBSERVATIONAL EFFECTS 

Observable effects consist of certain time variations 
of the signal during the nights, A^it). Variations have 
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FIG. 3: Effect of averaging over the production region of ^Be neutrinos in the Sun. The relative change of the electron 
neutrino flux for mantle crossing trajectories with rj = 0.58... 0.59 (left) and 1.5 ... 1.51 (right) without (dotted line) and with 
(dash-dotted line) averaging. 


quasiperiodic dependence on time, the amplitude of vari¬ 
ations also changes during the night, furthermore os¬ 
cillations with high frequency are modulated by lower 
frequencies. Sinusoidal dependence is modified by non¬ 
sphericity of the Earth, small structures of the profile, 
etc.. Parameters of these variations depend on position 
of a detector. 

In our estimations of sensitivities to different effects 
we will use for simplicity the spherical Earth. Devia¬ 
tions from sphericity should be taken into account for 
different detectors individually. As function of L the os¬ 
cillatory curve has regular dependence which is similar 
to the dependence on cos 77 for spherical case. Local de¬ 
viations from sphericity produce complicated effects and 
perturbations of the profile should be taken into account. 

In the case of spherical Earth and without small struc¬ 
tures we can use z = cos ry as variable and evaluate num¬ 
ber of events in small Az bins. Recall that different de¬ 
tectors cover the nadir angle ranges with minimal rjmin 
which depends on the latitude of the detector. 


Rate of events in the detectors 


The rate of — e scattering events in a detector can be 
written as 

R{t)=N, f dEFBeiE,t)a,iE)[PiE,t)+rNcii-PiE,t))], 

(44) 

where Ng is the total number of electrons in a fidu¬ 
cial volume, CTe is the cross-section of Ug — e scattering, 
r^c = o'Nc/<^e, and ctatc is the cross-section of the v^ — e 
scattering due to neutral currents. The last term in (44) 
accounts for the contribution of the and v^. neutrinos 
that are generated by oscillations. In the spherical Earth 
approximation dependence of the probability on time is 


via 77 (f) only: P[E,t) = P{E,r](t)) The rate (44) can be 
rewritten as 


Rjv(t) = NgF°^{t)dg f dEg{E) [P(e;, f)(l-r^c)+rivc], 

(45) 

where cr is the averaged (over the ^Be line) value of cross- 
section. 

Let us introduce the relative variations of the rate of 
events due to the Earth matter effect 


RN{t)-RD{t) 

= Ro{t) 


(46) 


Here Roit) is the rate which would be without Earth 
matter effect and it is given by Eq. (45) with P/v substi¬ 
tuted by Pd- Inserting (45) into Eq. (46) and assuming 
that 0 ^, and consequently, Pd do not change with en¬ 
ergy in the interval of energies of the ^Be neutrinos, we 
obtain 


Affit) = Ag{t)K, (47) 


where 


1 - rNC 

1 + r^c(l/PD - 1) 


(48) 


gives the correction due to contribution of and Vr ■ For 
rNC = 0.2 and Pd ~ 0.6 we obtain k = 0.7. Thus, the 
correction leads to damping of variations of signal. 

For estimations we will consider future scintillator (or 
scintillator uploaded) detector with fiducial mass Md = 
100 kton and 5 years exposure as illustrative values. By 
simple re-scaling one can find effects in specific proposed 
detectors such as LENA, JUNO or WBLS. 

To evaluate the total number of events in such a detec¬ 
tor we will use the rate estimated in (see also |25)): 
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1.5 • 10^ events per day in 48 kt. Normalizing to this 
number we obtain 

= 5.7 • 10^ f ^ f events. (49) 

\ lOOkton J \ 5years J 

These events can be analyzed in various ways to make 
the data sensitive to different quantities. 


Establishing the Earth matter effect 


For this it is enough to determine value Ajy averaged 
over rj. So, one should compare the total numbers of 
events detected during the nights and days during whole 
the exposure. (Possible distortion of Aj^{t) by perturba¬ 
tions of the density profile is not relevant here.) We take 
that approximately half of the total number of events 
(49) is detected during nights and another half during 
days: 


Nn^Nd = 2.85 • 10^ events. (50) 

Then difference of the numbers of night and day events 
due to oscillations in the Earth equals 

No — Nm « —NoAjq = —NonAe = 1.7 • 10"^ events. 

(51) 

The statistical error of measurements of No—Nn is ctat = 
^/2Nn = 7.5 • 10^. So, after 5 years 


{ND-NN)^2.3aN. (52) 

That is, the difference of the night and day signals can 
be established at 2.3 a level. This is in agreement with 
estimation made in the introduction. Notice that various 
systematic uncertainties cancel in the relative variations. 


Variations of the signal during nights 


One can measure oscillatory variations of the ^Be neu¬ 
trino signal during nights (Fig. detecting events in 
short time intervals. This method may have less system- 
atics than the first one described in the previous sub¬ 
section. However here new systematics may appear due 
to unaccounted effects of small structures of the Earth 
(see below). Period of time variations can be estimated 
for spherically symmetric profile in the following way. 


The period in the nadir angle scale is given in (40), 
rjT = 2.24 • 10“^/sin 77 . Then the period in time equals 


tr = 




drj 

2Re sin 77 dt 


-1 


(53) 


where the speed of change of the rj with time depends 
on value of rj. We find that the average period equals 
about 1 min. So, to measure the oscillatory curve one 


needs to take time intervals smaller than 15 sec. Number 
of events expected in such an interval will be about 3 - 
5. Therefore summation of signals from time intervals in 
which the Earth matter effect is the same during several 
years is needed. 

In the approximation of spherically symmetric Earth 
the length of trajectory and the density profile are fixed 
uniquely by the nadir angle 77 , and the Earth matter effect 
is quasi-periodic function of cost? (since L = 2Re cosrj). 


Then from (34) we obtain the period in cos 77 


(cos 77) T = 


L 


2Re{1 — cos 20 i 2 cf^e )' 


(54) 


Due to increase of 6 ( 77 ) with decrease of 77 the period 
will slightly increase for deeper trajectories. The relative 
change of period is of the order 10“^. It would be exactly 
periodic function in the case of constant density. 

So, the method consists of splitting the whole cos 77 in¬ 
terval for a given detector into small intervals A(cos 77 ) < 
( 00377 ) 7 ’, identification the corresponding time intervals 
during nights. At = At(cos 77 , A(cos 77 )), and accumula¬ 
tion of events in these small intervals during several years. 

However, even with 100 kton detector it is not possi¬ 
ble to measure whole the oscillatory curve of Eig. [^with 
appreciable statistical significance. Therefore one should 
also sum up the signal over all periods of Fig. Intro¬ 
ducing corrections due to change of period with 77 one 
can combine all the events in a single period using the 


effective phase (34) as variable. This means that for each 


event or several events detected during small enough time 
intervals At = (10 — 15) sec, one finds ^^(t), and then 
collects events for the intervals ^rn(t) + 27rk. 

Following this we divide all the events detected during 
the nights in to two groups: events detected in the first 
half of period, Ni , and events detected in the second half 
of period, N 2 . Then difference of the events equals 


iVi - 7 V 2 = (D}Nn « (D}No, (55) 

where 

(D) = (56) 

TT 

Here A^°'^ and maximal and minimal values of 

An (averaged over different periods) and factor 2/tt re¬ 
flects decrease of the difference due to integration over 
half periods as compared to the total depth. Thus, (D) 
is the relative variation of number of events integrated 
over half a period. From Fig. we find that for the 
mantle trajectories — A^™) « 5.2 • 10“^. Taking 

Nm ~ 2.8 • 10^, we obtain Ni — N 2 = 0.9 • lO'^. The sta¬ 
tistical error of measurements of the difference {Ni — N 2 ) 
is VNN = 5.3 • 10^. Correspondingly, variations during 
the night can be established at 1.8cr level. 

Let us estimate effect from the core crossing trajecto¬ 
ries. We assume that fraction of events collected from 
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FIG. 4: The relative change of the electron neutrino flux 
for the mantle crossing trajectories as the function of rj for 
two different values of width of the ^Be line which correspond 
to two different temperatures in the center of the Sun: Solid 
line for sun central temperature Tq = 15.55 x 10® (solid line) 
7.77 X 10® degree (dashed line). 


these trajectories is about 10 %, that is, Ncore = 2.8 • 10 ®. 
The depth of modulation is larger: = 

2 • 10“®. The expected difference of events in the first and 
second halves of the period equals Ni — N 2 = 2.5 • 10^. 
The statistical errors is cr = '/Wjj = 1.7 • 10^. Therefore 
Ni — N 2 = 1 . 5 ( 7 . Summing up the significances in the 
core and mantle we obtain that variation of the signal 
can be established at 2.3 (t level. These estimations have 
been performed for the ideal periodic oscillatory depen¬ 
dence of Ajq on COST?. As we discussed in Sec. 3, the 
presence of small structures at the surface of the Earth 
and in the crust distorts the sinusoidal dependence, in 
particular it can shift significantly maxima and minima 
in the cos 77 scale. This leads to uncertainties in determi¬ 
nation of borders of periods and intervals in which events 
should be summed up. The uncertainties can be reduced 
if a profile of the Earth along trajectory of each event is 
known and therefore for each event the oscillation phase 
can be determined. Actually, it will be enough to know 
the length of trajectory, and the rest can be accounted as 
systematic error. As a result, significance of establishing 
of variations may be somehow lower. 


Determination of the line width 

The depth of oscillations in the Earth decreases with 
increase of the width of the line which in turn is propor¬ 
tional to the central temperature in the Sun: Tee oc Tc- 
In Fig. I^we show the oscillatory picture for two different 
widths of the line which correspond to two different cen¬ 
tral temperatures of the Sun. As follows from Fig. the 
depth of oscillations increases by factor 2 with decrease 
of the temperature by the same factor. 


According to Fig. [^for shallow trajectories, 77 > 1.4, 
the averaging effect is negligible and the depth of oscil¬ 
lation is the same as for monochromatic spectrum. The 
averaging effect strengthens and depth decreases with de¬ 
crease of 77 . For larger width Tee, the change of depth 
is faster. The depth of oscillations is larger for the core¬ 
crossing trajectories and it changes with 77 weakly. There¬ 
fore methods of determination of the width can be based 
on determination of the depth of variations Aj^ and its 
dependence of the depth on 77 . One can 

1 . measure the average depth of variations and com¬ 
pare it with expectations; 

2 . compare the average depths for deep (e.g. 77 = 

0.58 — 1.20) and outer (77 = 1.20 — 1.57) trajectories in 
the mantle; 

3. compare of the average depths in the core and in 
the mantle. 

Let us consider the first method. According to Fig. 
the depth {D) oc ^ and 11 ( 77 ) changes linearly with 77 (see 
Fig. 0 . Therefore we can write for the averaged depth 
for all mantle trajectories 

0.19^) • 10-®, (57) 


{D) = 


0.63 


where is the standard value of temperature in the 
center of the Sun. The first term in the brackets cor¬ 
responds to the half depth of oscillations at the surface 
(trajectories with large 77 ), the second one corresponds to 
the deepest trajectories in the mantle for Tc = T®*. For 
Tc = Tc* we obtain from (57) {D) = 0.51 • 10“®, whereas 
for Tc = 0.5T®L (D) = 0.64 • 10“®. So, decrease of the 
Tc by factor 2 leads to increase of the average depth by 
20%. As we have found in the previous subsection the 
deviation of the average depth from 0 can be established 
with significance about 2.5 (t, so its 20% change would 
correspond to about 0.5a. 

Consider the second method. According to Fig. [l] the 
average depth in the outer part of mantle (77 = 1.20 — 
1.57), equals (D) = 0.64-10“®, whereas in the inner part 
{rj = 0.58 - 1.20): (D) = 0.33 • lO"®. We take that the 
number of events in the inner and outer intervals of 77 
are the same and equal Nj^/2, so that variations in the 
outer and inner ranges of 77 equal 6.7 • 10^ and 3.4 • 10^ 
correspondingly. The statistical error is 3.9T0^, therefore 
variations can be established at 1.7a (outer) and 0.9a 
(inner). The ratio of the depths (which is sensitive to 
temperature) equals 


{D)iT, 

{D)ou 


= 0.51 ±0.66, 


where the Icr statistical errors is indicated. So, even de¬ 
crease of the ratio by factor 2 will correspond to 0.4a. 

Sensitivity is low but this is independent measurement 
of the width. Clearly, factor of 2 uncertainty in central 
temperature is excluded, e.g., by measurements of the 
boron neutrino flux. 
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The oscillatory pattern of Fig. depends on precise 
value of Am 2 i- Uncertainty in Arri^i, S^Am^i), would 


not influence the pattern if 


djAmli) 


< 


0.3%. 


(58) 


The present accuracy of determination of Am^i is 2.5% 
and future experiments (e.g. JUNO-reactors) will not 
be able to reach 0.1% accuracy. So, Amfi should be 
extracted from the ^Be-neutrino studies simultaneously 
with other measurements. 

The mass splitting Am^i determines the depth, the 
average value of probability and the length of oscilla¬ 
tions. Correspondingly, there are several different meth¬ 
ods to measure Arri^i- The depth of variations of 
and its average value according to Eqs. ( [I7| ) and (21) 
are inversely proportional to the splitting: D oc. 1 / Am^, 
Af. oc XjArn^i- As we have established in the previous 
subsection, A^ and D can be distinguished from 0 at 
about 2 . 3(7 level. Consequently, these measurements will 
determine Am^i with about 1/2.3 « 40% accuracy at Ict 
level. 

Since there are many (up to 400, according to Fig. 
periods of oscillations, dependence of the oscillation 
period on Ato^i can provide much better determination 
of Am 2 i- Indeed, the error 6{Amlj^) ~ (l/300)Am2i 
would lead to substantial change of the oscillatory pic¬ 
ture. Clearly knowledge of exact values of L for the de¬ 
tected events is crucial. 

The analysis can be performed using, e.g., the Lomb- 
Scargle m, EH] method which is valid for uneven time 
intervals. Let 


Nf, = NUl + A^M) 


(59) 


be the observed number of events in the j— bin with 
the average nadir angle ijj. This corresponds to the true 
value of Am 2 i = Am^™® The first term in (59) can be 


considered as fluctuating noise, whereas the second one 
as periodic signal to be extracted. For constant density 


we can use expression ( 21 ), so that 


obtained by varying Pls over Amlf it corresponds to 
maximum of Pls- 

For illustration we can perform simplified computa¬ 
tions using constant density profile and spherically sym¬ 
metric Earth. In this case the N^- is exactly periodic 
function of cos 77 . 

The Lomb-Scargle periodogram is based on the dis¬ 
crete Fourier transform, but in the case of many bins of 


even size we can substitute summation in eq. (61) by 


integration over x oc cos 77 from 0 to x^ax (here x is the 
length of the trajectory). So, 


Pls = 
where 


[Xrr^ax I dx[Njv(x) - Nq] COS A^^xf^ 
XiXmax J" dx COS^ Ajji X 


1 


Nn{x) « A£,(x) ( 1 -b -cjge/cos A‘ 


The integration in (62) gives 

(chef? 


Pls = Nd ' 


sin A$ 


A$ 


where 


A* — - AP*)r 

JXr, 


(62) 


(63) 


(64) 


(65) 


Using expression for A^ in vacuum we find the follow¬ 
ing relation for S{Am 2 i) = Attt,^™® ^ — Am^i 


S{Am2i) A$ 


Amj^ 


AtT Xn 


( 66 ) 


For Ad) — TT12, which corresponds to half of the height 
of the peak in (64), Eq. (| 66 | gives 


SjAml^) 

Aml^ 


= 4-10 


-4 


(67) 


These estimations do not take into account effect of av¬ 
eraging over the ^Be-neutrino spectrum. So, realistic ac¬ 
curacy will be above 0 . 001 . 


1 

2-13- 


fV^«7V|,(l+ic?3e/cosAr^L,). 


(60) 


Tomography of the Earth 


The method require that the mean value of the back¬ 
ground is zero. So, we need to subtract the average 
day signal Njy from Nj^. Then the Lomb-Scargle pe¬ 
riodogram is defined as 


Pls = 


1 - Np] cos AU^ai,)^ 

n X), cos^ 


(61) 


where n is the total number of bins and A^* corresponds 
to the fit value of Atti^i. The true value of Am^i can be 


As we discussed in Sec. 3, the oscillatory pattern en¬ 
codes information on the shape and density profile of the 
Earth. So, in principle one can perform tomography of 
the Earth with spatial precision comparable with the os¬ 
cillation length. Let us summarize some dependences: 

- Sudden change of the depth of variations of Ajv at 
77 ~ 0.58 marks trajectories which start to cross the core. 
The corresponding 77 gives the position of the density 
jump between the mantle and the core, the size of change 
of the depth of oscillations reflects the size of the jump. 
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- Detailed oscillation pattern in the core region is sen¬ 
sitive to the parameters of the inner core. 

- Small density jumps in the mantle trajectory range 
and modulations of the high frequency oscillations en¬ 
code information about the mantle structure: position 
of borders between different regions (inner mantle, outer 
mantle, etc.) and sizes of density jumps at their borders. 

- The oscillatory picture is sensitive to the shape of the 
Earth, in particular, to its deviation from sphericity. 

- Distortion of the periodic oscillatory curve is sensitive 
to local structures (mountains, seas, oceans, oil layers, 
etc.) [TT] . 

Detectors with 0.1% accuracy of measurements will see 
(resolve) core and mantle. Other features will have lower 
significance and appear as systematic errors in these de¬ 
tectors. 


5. SEARCHES FOR STERILE NEUTRINOS 


The parameter e, and consequently, the Earth mat¬ 
ter effect are inversely proportional to the mass splitting 
Therefore with decrease of Amfi (for the same 
value of mixing) the variation Ajv increases. In this con¬ 
nection let us consider sterile neutrinos with very small 
mass splitting, Am^Q = ttiq — ml <C Am^^- Here TOq 
is the mass of new state. The matter parameter for 
the Ve — Vs system, e,,, is determined by the potential 
Vs ~ 0.5Ve. Then the resonance, Cs = 1, is achieved for 

Amlo = = 0-9 ' 10”^ eV^. (68) 

(for the mantle densities). In the core the resonance en¬ 
hancement will be at about 2 time bigger splitting. 

For definiteness we will consider mixing of sterile neu¬ 
trino in the mass state vi characterized by the angle 9s. 
In this case we obtain for the difference of the night and 
day signals 

AP « |C/eiP(|HoiP - |Hoon(Pfe - (69) 


where Aqi and Aqq are the amplitudes of transitions be¬ 
tween the matter eigenstates: i/Qm vi^ = vi and 
vom —t t'om = vq when neutrinos propagate from the 
center to the surface of the Sun. In (69) we have taken 
into account that in the Sun, for small mass splitting the 
neutrinos are produced far above vi — vq resonance and 
in initial state t'om ~ vi which gives the factor |C/eip in 
Eq. (69). Eq. (69) is analogy of expression (10). 

Even for splitting ( 68 ) the coherence will be lost due 
to separation of the wave packets of the mass states on 
the way to the Earth. So, as in the active neutrino case, 
independent fluxes of vi and vq will arrive at the Earth 
surface. 

For oscillations in the Earth we consider the 
2 i/—system vi — vq and denote by Pfg the probability 


of vi —>■ Vf, transition. For estimations we use the con¬ 
stant density profile for the trajectories inside the Earth 
with average value of density. This is justified since due 
to smallness of Ato^i Hie oscillation length is much larger 
than the size of the Earth. In this case 

- |C/ei|" = Sin" 2C sin" Qa^L^ , (70) 


For relatively large mixing angle the adiabatic transition 
vom —t Vq inside the Sun would lead to substantial ad¬ 
ditional suppression of the Vf, ^Be-neutrino flux. Theo¬ 
retical accuracy of the flux is about 1.4%. Experimental 
accuracy (presently 5%) will be much better. So, the 
survival probability should be larger than 0.98, that is, 
the adiabaticity should be strongly broken with |Hoi|" > 
0.98. Using results for active neutrinos (see e.g. [22]) 
we estimate that this can be achieved if sin" 20s < O.OI 
for splitting in ( 68 ). Taking |Hoi|" — |Hoo|" « 0.96, 


|UeiP ~ 0.67 and Pd ~ 5/9, we obtain from Eqs. (69) 


and (70) 


Ajv ~ —l.I5es sin" 20 ™ sin" ( -AmL 


(71) 


Let us consider maximal allowed value sin" 20s = 0.01. 
In this case the width of the MSW resonance, AE = 
2tan20sE « 0.2E = 190 kev, is much larger than the 
width of the line. Changing Am"Q within (20 — 40)% 
would put the line at different points of the MSW reso¬ 
nance peak. 

In the resonance we have Cs ~ sin" 29 sm ~ 1 and there¬ 
fore 


-1.15sin" l-AmL{r]). 


(72) 


For the resonance value of Am^i ( 68 ) the oscillation 
length in vacuum equals l^, = 2.5 • 10"^ km and in mat¬ 
ter (in resonance) Im = ?i//sin 20s = 2.5 • 10® km, i.e. 
much larger than the diameter of the Earth. For the 
deepest trajectory in the mantle, L « 10"^ km, the phase 


equals 0.126 rad., and according to Eq. (72) Ag = 0.017. 


This is an order of magnitude larger than the effect for 
active neutrinos and of the order of the present accu¬ 
racy. For cos 77 = 0.4 (middle trajectory in the mantle) 
we hnd Ag = 0.0045. The effect monotonously increases 
as oc l/(cos 77 )", reaching maximum for the deepest tra¬ 
jectory (middle of the night) and then it decreases down 
to zero. 

For the core crossing trajectories the effect can be more 
complicated. The MSW resonance in the core leads to 
appearance of another peak, and the ^Be neutrino line 
can be in one peak or another. Also the interplay of the 
effects in the mantle and the core may occur. 

If Amfg is more than 40% larger than the resonance 
value, the values of oscillation parameters become close 
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to vacuum values: Im ~ Iv, sin^ ^Osm ~ 0.01, so that 

Cl ~\ 

Ae « -0.015es——sin^ -/\^L{r]). (73) 

In this case still £« = 0(1), and the oscillatory factor can 
be of the order 1. Thus, the effect for the deepest trajec¬ 
tories in the mantle may reach 1%. Since the oscillation 
length becomes comparable with the size of the Earth, 
the time dependence of the effect is more complicated, 
e.g. with two maxima symmetrically shifted from the 
middle of the night. 

With further increase of Am^Q (outside the resonance 
region) the size of the effect decreases as l/Am^g, and 
time dependence will acquire an oscillatory form with 
increasing number of periods. 

Let us make similar estimation for sin^ 26s = 0.001. In 
resonance the oscillation length equals Im = 8.5 • 10® km 
(about 3 times bigger than in the previous case). Corre¬ 
spondingly, the phases will be 3 times smaller, and the 
oscillatory factor will be an order of magnitude smaller 
As a result, for the deepest mantle trajectory Ag ~ 10“®. 
Outside the resonance peak we obtain < 0.001. Thus, 
the effect decreases as sin^ 20s everywhere. 

6 . CONCLUSIONS 

We explored in detail effects of propagation of the solar 
^Be neutrinos in the matter of the Earth. We estimated a 
possibility to detect these effects with future large scintil¬ 
lator (or scintillator uploaded) detectors which will have 
0 .1% accuracy of measurements. 

The main features of the propagation are determined 
by low energy of the ^Be neutrinos and their narrow en¬ 
ergy spectrum. Oscillations in the Earth are pure matter 
effect and they are related to transitions between the 
mass eigenstates. 

Physics of oscillations is determined by two accidental 
coincidences. Due to low energies, the oscillation length, 
« 30 km, is of the order of small structures of the Earth 
profile and non- sphericity of the Earth. The width of 
the spectrum is comparable with the period of the oscil¬ 
latory curve in the energy scale. So that depending on 
length of trajectory (nadir angle) one should observe dif¬ 
ferent degree of averaging. In configuration space this is 
equivalent to partial loss of coherence due to shift of the 
wave packets of different eigenstates of the Hamiltonian 
in the course of propagation inside the Earth. The size 
of the packets due to spread on the way from the Sun is 
several orders of magnitude larger than the shift. 

The main observable is the oscillatory variation 
(mainly suppression) of the signal in time (with nadir an¬ 
gle) during the night. The depth of oscillations changes 
with time since the averaging of oscillations becomes 
stronger with increase of the length of trajectory. The 
depth is the largest for shallow trajectories and it is the 


smallest for the deepest mantle trajectories. It sharply 
increases when trajectory crosses the core. Small density 
jumps in the mantle and the core produce modulations 
of the oscillatory curve. 

We find that for illustrative configuration of experi¬ 
ment (100 kt, 5 years of exposure) the Earth matter ef¬ 
fect can be established at 2.3cr level; the width of the 
^Be-neutrino line can be determined with factor of 2 ac¬ 
curacy at 0.5 (t; Am^i can be measured with accuracy 
0.1%. The presence of the core of the Earth and its bor¬ 
der can be seen by the detector (close to equator) at 2cr. 
The sensitivity can be enhanced if one uses larger expo¬ 
sure time (e.g. 10 years) or larger fiducial volume of the 
detector. 

Determination of other characteristics: size of region 
in the Sun where the Be neutrinos are produced, detailed 
tomography of the Earth (small scale structures, layers in 
the mantle and the core) will require further substantial 
increase of the detector size. 

At this level of accuracy it is not possible to measure 
whole oscillatory curve which contains interesting infor¬ 
mation about structure of the Earth. That would open 
up a possibility to use the Sun as the scanner of the Earth 
to perform tomography. However, estimations show that 
the required sensitivity is not by several orders of mag¬ 
nitude higher. 

One can perform searches of sterile neutrinos with mix¬ 
ing sin^ 26s = 10“® — 10“^ in wide range of Ato^q > 10“® 
eV^, especially in the resonance region of mass splitting 
around 10“^ eV^. 
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